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PREFACE 


The technical results presented here represent a ten week 
effort by its authors during the Summer Institute for Bio- 
medical Research sponsored by the Technology Utilization 
Office at the Goddard Space Flight Center. Their chal- 
lenge wus to apply NASA developed technology toward the 
solution of this particular problem and to demonstrate its 
usefulness to other problems in medical diagnostic moni- 
toring instrumentation. 

This report has been puulished and made available for 
general use so that others in both the technical and med- 
ical communities might benefit from the work of these 
individuals. 



Summer Institute for Biomedical Research 
Technology Utilization Office 
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June 29, 1970 


PROPOSED PROJECT 


After reviewing the previous work, either design a new circuit to produce heart 
sound envelopes, or improve the existing one. Design requirements are to pro- 
vide frequency information to some extent, calibrate the signal for true intensity 
and produce a portable prototype, preferably of very small size. The circuit 
which meets these requirements should thereby make it feasible to identify either 
by physician or computer the interpretations listed on page of this notebook. 

Routine clinical evaluation of the phonocardiographic amplitude/time signal is 
subject to interpretation of the complex waveform which carries a full continuum 
of frequency information. The true intensity and frequency relations are very 
difficult to perceive in this situation. A circuit which would produce an intensity 
(power) signal from the amplitude (voltage) signal would provide information 
better related to what the physician heard, i.e. , loudness of sounds and their 
temporal relations. Frequency related information, analogous to pitch of audible 
heart sounds would need to be retained to some extent, but not to the equi^ alent 
of spectral analysis methods. 


Discussion with Dave Winer from George Washington University on the sound 
envelope project. Topic: Preprocessor for Health Soimds Output of Preproc- 
essor: instantaneous value of intensity with respect to time. 


The output of the preprocessor must have the following characteristics: 

1. Gives the physician a picture of what he hears 

2. Gives the computer a waveform that can be easily stored and analyzed. 

The output of the preprocessor can be shifted in time from the real soimds— 
although it is desirable not to have a time shift. 


In the final design the preprocessor should be small enough to connect to the 
microphone or to a recording unit. 


July 1, 1970 


Notes from: Winer, et al. ; Heart Sound Analysis: A Three Dimensional Ap- 
proach; The American Journal of Cardiology, Vol. 16, October, 1965 

"It is difficult to relate audible intensity of complex sounds to instantaneous 
values of amplitude as portrayed in the osuillogr^hic form. " 
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”Visible lluctuation in amplitude may not in many instances eoia espond to audi- 
ble sensations because the eaj- resjjonds to root- mean-square pressui'e, which 
means that sound must persist Tor some duration in oi’der to be perceptible. 
Some transients visil)le on the oscilloj>raph, thcreiore, may be oi insuiiicient 
duration for audible recognition. " 

Notes from: Winer, David; Notebook: Phonocardiograin Project, Deceml)er 1, 
1964. 

”A simplified wavefoi’m c:m be obtained which will include the same inloianation 
which is presently used and needed. But it slumld have these adv antages: 

1. easier to program 

2. easier to interpret visually 

3. slower sampling rate 

a. compres‘-:ed feed time to com[)uter 

b. data phone tnmsmission" 


Method 


"Before the PCG signal is digitized it eiui be full wave I'ectiiied. The ne.xl stage 
in the circuit would be introduction of a suitable time constant to prev ent return 
to l)aseline (this constant can be adjusted for the desired resolution loi’ timing 
of split sounds, etc.) This is essentially the technique descri'oed liv:” Uushmer, 
et al. ; "Sonvelographic Recording of Murmurs During .\cute Myocarditis:" 
Americim Heart Jounial, 48:835, 1954. 

"The result is still an amplitude curve, wducli does not divulge intensity. 11 one 
wished for an intensity /time curve, this could l)e accomplished by a circuit sim- 
ilar to that in a vacuum tube voltmeter which has output calibrated in decibels. 
This output signal could be put on tape at the data acquisition system or could be 
generated at the computer center immediately prior to digitizing.” 

Note added; (February 19, 1965) 

"Si>eech data are carried lai’gely by the varying .shape of the |x>vver density spec- 
tniin and not— as many wrongly believe— in the sound pressure vs time plot seen 
on an oscillosco|x'." 


S. K. Gerber and E. J. Strausman ("speech scientists," 

Communication Division Hughes A/C Co. , 5440 W. Century, 

Los Angeles, California) in Digital coding for wire commu- 
nication, Space/Aeronautics v. 40, No. 5, October, 1963, 
pp. 95-96. 

(February 25, 1965) "A baseline was determined by finding the MODE of all val- 
ues of amplitude. (The mode is that value which appears most frequently.) We 
have found that this makes an excellent baseline - Don Sherman’s contribution. " 

(July 29, 1966) "If an analog network were designed to create this type of output 
(visual records of the way heart-sounds sound to the ear) from the microphone 
amplitude signal, it could be handled much the same as other low-frequency sig- 
nals, such as ECG, PTG, EEG, etc. We could expect to find the following 
benefits: 

• 500 digitizing rate vs 3000 

• FM channel recording 

• Easy computer pattern recognition program 

• Telemetry over standard circuits 

• Lower pulse code modulation rate if NASA would like to send heart- 
sounds from space 

• Less computer care" 

(November 31, 1968) Referring to Mark Wilber's circuit: 

"This technique can be used to 'preprocess' heart sounds before recording, al- 
lowing all data acquisition, telemetry, preprocessing and processing the com- 
puter using the standard MSDL ECG hardware." 

(December 4, 1968) Comment about computer analysis: 

"Note that this is a digital computer approach to the intensity curves of pp. 6, 

7, and 8 of this book. The analog preprocessor saves MUCH computer time." 

Perry et al. ; Computer Analysis of the Phonocardiogram; reprinted from Engi- 
neering in t he Practice of Medicine ; The Williams and Wilkins Company, 1967. 

"Rectification and smoothing techniques may aid in more precise identification 
of the heart sounds in the presence of artifacts or murmurs. " 


July 2, 1970 


Siebert; from notes given in course titled "Signals and Systems" given at MIT in 
Spring, 1969. Notes were written in 1967. 

Square -law and other non-linear-res istive devices 

An ideal square-law device is described by the formula 

Y(+) = A (+) (2) 

(where A is a constant) and is often used to represent approximately the input- 
output characteristics of full-wave rectifier circuits such as shown in Figure 1. 



Figure 1. Two Examples of Full-Wove Rectifier Circuits 


Other common devices which might be represented by (2) are the temperature of 
the heater element in a thermocouple-type ammeter as a function of the current, 
and the output of photomultiplier illuminated by a laser. The actual input -output 
characteristics of such devices might more accurately be described as in Figure 
2. The extent to which such a graph will be adequately represented by (2) de- 
pends, of course, on the particular diodes employed, on the amplitude of the 
signals and on the precision required in the representation. 



e, (+) 


Figure 2 


Idea # 1 
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time constant of envelope detector = RC 

The voltage across the capacitor may have an undesired ripple frequency which 
may be filtered out by a low-pass filter. 

Diodes used for full wave rectifier 


I (mal 



V (volts) 


rectifier diodes and D1 : DR435 diodes 
C = 20ufd, electrolytic capacitor 
R = rEtnges from 0 to 5 M ohms 

changing the value of R the time constant varies from 0 to 20 sec. 

Problem with the four diode full wave rectifier: 

It is desirable to have both sides of diode D2 above grounded. With both sides 
of this diode grounded there will not be full wave rectification making it neces- 
sary to find another full wave rectifier. 
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lion Notes 



August 27, 1970 





Full Wave Rc'ctilier Used: 


2K 



This lull wave rcctificT is JR iiH of the frequencies that are needed to rec- 
tify heart sounds. It has a low output impedance anti there is no ^roundinn prob- 
lem which makes it a ^ood match for the next sta^e of the preprocessor. 


Aujjust 2H, 1970 


Knveloix* Detector: 



This enveloix' tietector simply smotjths out the rectified siijnitl. H does not lilter 
out any part of the signal that Is needed retaining the information from all the 
frequencies contained in a heart sound. 

An RC time constant is chosen so that the necessary diagnostic information is 
retained while at the same time the waveform is smoothed. Time const:uit 
used - RC = 0.02 to 0.03 seconds. 


On the iKJsitivc cycle of the input sinn;il, the capacitor C (lOul) ehart;es up to the 
pc*ak voltJige of the input sifpial. As the input signal falls Ix'low the |x*ak value, 
the diode is cut off Ix'causc the capacitor voltJige (which is very nearly the |X‘ak 
voltJige) is greater than the input signal voltage, thus causing the diode to ojxn. 
The capacitor then discharges through the resistor R. During the next jxisitive 
cycle, near the peak of the input signal, the input signal Ixconies greater than 
the capacitor voltage 'Uid the diode conducts. The capacitor again charges to the 
peak value of this ■new cycle. During the cutoff ix*riod the capacitor will dis- 
charge completely with no new input. 

IjOw Pass Filter 

"Frequency selective networks for use in the frequency range below lOOkllz have 
always been a problem. In this area of operation the inductors and capacitors 
required are large, Ixjth in value and physical size. Also, at these fi’equencies 
inductors and capacitors Ixjcome quite lossy and the circuit Q's tx*gin to suffer. 

"The answer to this problem is to exchange the large inductor ;u»<l capacitor for 
a large block of gain and use well knmvn feedback principles to achievt' selectiv- 
ity with R-C active filters. Previously, to achieve a high degree oi accuracy 
and circuit stability, a large numKr of active components was n quired in a 
fairly sophisticated circuit. Consequently, the design time and num>xr of active 
components required made the use of active filters quite exixmsive. 

"The solution to this problem came with the advent of integrated circuits which 
allowed transistors to lx "less expensive" thim resistors. Now, excellent gain 
blocks cun lx* fabricated at fairly reasonable costs. .Vnd as technology improves, 
the ixM-formance will continue to improve and the costs will continue to decline, 
making the use of active filters very economical." 


from Motorola Application Notes 

.\f er rectifying the heart sounds and putting the rectified wave into the envelope 
detector, we put the output of the envelop*’ detector into a Rockland filter to de- 
termine what cutoff frequency was needed for a smooth envelop*’. Thi.s filter is 
filtering the envelope, not the heart sounds . Desired cutoff fre*ii'viicy - 
cycles 
sec. 


K 


Second Order Low Pass Active Filler 



where 


d 


2 + R|C: 
1 2 



1 - R, C, 

2 = R: C: 


want 


n = 35 — so Wn = 2 


sec. 


rad. 

sec. 


pick R] = 120 K and C] - 0.05 

1 - R|C| =(1.2xl0') (0.05x10'^) - 6x10"^ sec. 


Wn- - 2 = 4.84 (— Y 2x lO^* 

12 (6x 10 ) \sec./ 


2 


1 

(G X 10*^) (4.84 X 10“ ) 


3.4 X 10"^ sec. 


pick 


R> - 13K 


r 2 — R iCi 


C: 


2 


3.4 X 10'^ 
1.3 X 10^ 


0.2G 
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Low pass active' filler used for preprocessor (from calculations on previous page) 



Operational amplifiers use: 
Model 4009 Medium 
Voltage D.C. Operation 
Amplifier from: 

OPAMP lABS 

172 S. Alta Vista Blvd. 

Los Angeles, Calif. 9003G 
(213) 934-356(3 


Cutoff frequency determined in the lab 


33 

sec. 


Latest full-wave rectifier used: 


2K 



Operational amplifiers used: Burr Brown 


August 20, 1970 


Note from Innocent Murmurs 


.Murmurs - 

The intensity of the murmur is related to the velocitj' of blood flo\>^ and according 
to one theory intensity varies as the fourth power of velocity of flow. 

Velocity is dependent on factors like volume of shunt, cardiac output, etc. So, if 
a small shunt and low' velocity of blood exists the murmur will lie of low intensity. 
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August 28, 1970 


AM Approach for "Demodulation" of Heart Sounds 


From the meetings I had with Mr. Mark Wilbur the following points l)ecame 
clear about this method and the way it had been used. 

1. Assume heart sound is an amplitude modulated signal and try to demod- 
ulate it. 

2. Pick a suitable carrier frequency to do the demodulation. 

3. Pick a suitable filter (cutoff frequency and voll. off) for blocking re- 
dundant information. 

4. Devise a method to multiply the carrier into the AM signal. 

The carrier frequency was picked at 20 Hz in order to demodulate the signal by 
sending the lowest end of the audible range of the spectrum. Any frequency 
lower than 20 Hz cimnot be heard, therefore, we do not (it was assumed) need 
that part of the frequencies of the signal in the envelope. 

The filter was picked as low pass filter with the cut off point at 17 Hz with a roll | 

off of 24db/octave. This choice was made by looking at the output of the filter \ 

and changing the cutoff frequency and the roll off until a reasonable envelope was i 

achieved. I 

The device to multiply the signal with the carrier signal was chosen to be a I 

diode . I 

The simplified circuit which was used to get the result was the following. * 



The assumption was tliat the diode in the above circuit could behave as a multi- 

V' 

plier. Current in the forward direction of the diode If = V ■*" 777 Tf, 

M • O • • • • 

using the exponential characteristic of a diode where V is the voltage across it. 
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I( - (\'j + cos 27rfot) 
If — Vj + cos 27rl'ot 
If ^ V, + cos 27rf()t 


(Vj ^ cos 27rfot)- 
2 ! 

+ Xl + cos - 27Tt\,t ^ 
2 2 

Vj" 1 -t cos 4jrfot 


(Vj + cos 27rfj,t)’ 

3 ] ^ 

cos 27Tl ot Vj 
+ V, cos 27rf^,t 


This current goes through R | and converts itself to a voltage value to go into the 
high iniiX'dance Op Amp to eventually go through the filter. The filter having a 
cut off frequency at 17 Ilz blocks the cos(2v t) obviously, blocks cos (4n t) 

V - V ■ 4 

blocks most part of (because is formed of higher frequencies because 


of the fact that most of Vj spectrum is alxive 20 Hz, so most of 



will be 


above 20 Hz, ) Most of the Vj term will be blocked too, so we will end up with 
the Vj. 


The results which were obtained from this method were OK because the output 
of the filter looked like the envelojje of the ineoping signal with approximation. 


The high frequency murmurs were not being detected and the low intensity por- 
tion of the signal was not detected well. The important point was that the pa- 
tient’s murmurs were usually of higher frequencies and sometimes of low inten- 
sity. So, the diagnostic value of the envelope was almost very small, if any. 


In order to get better results from this method we replaced the diode with a real 
multiplier. The multiplier vms a Fast Quarter-Square Multiplier/Divider. It 
was a definite improvement over a diode. The diode was picked on the assump- 
tion that we could ignore the higher terms in the Taylor Series of the exponential 
curve of the diode. This is usually not a valid assumption with the kind of sig- 
nal we have here. The voltage goes high and low and we know when V is not very 
small. The assumption of ignoring the higher terms is invalid. Furthermore, 
in order to get better results we decided to improve the choice of carrier fre- 
quency and also filtering. 


The Rockland Filter was used as a frequency analyzer to determine a spectrum 
knowledge (intensity vs time) of the heart sound. The bandwidth of 20 Hz was 
picked and the intensity of the sound was measured in this bandwidth as it was 
swept across all frequencies involved (0-1000 Hz). Many curves were obtained 
for different heart sounds. The conclusion was most of the spectrums were very 
different from each other, but the fact was that most of them had high intensity 
components (usually belonging to first and second sounds) Ijetween 40 Hz and 
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100 Hz. This knowledge determined the choice of carrier frequency for us. It 
was picked to Ije 40 Hz. This would shift the "important” part of the spectrum 
down to D.C. and around it. The filter cutoff was at 25 Hz with a roll off of 
24db/octave. The result was a signal similar to the envelope of the sound, but 
again very approximate and poor with respect to higher frequency murmurs. 

One important mistake with the previous circuit is that because of the diode haif 
of the sound signal is completely ignored. In any case, the second attempt with 
the AM approach failed too. This was lx?cause the spectrum of the sounds we 
achieved showed that no one frequency can be picked to demodulate the signal 
with, for the simple reasons that the complex signal received from the heai-t is 
not an amplitude modulated signal and it is not even close to something of that 
kind. The spectrums were very much distributed and there were no definite 
peaks to suggest choice of several carrier frequencies for "demodulation." 

This idea to pick different parts of the spectrum, shift them to DC and low fre- 
quencies and then add them to each other, fails to be a good one for several 
reasons. As mentioned above there is difficulty with choice of carrier frequen- 
cies. Then there is the problem of the negative frequency components of the 
Fourier Transform of the signal which gets shifted to DC and low frequencies 
each time. 


We multiply the signal with cosine of some frequency of time. 


INTENSITY 
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In (a) we have spectrum of a typical heart sound which is shifted down in (b) and 
the resultant spectrum is shown in(c). The overlapping of the negative frequency 
components introduces distortion of the results. This error could be very much 
amplified when several carrier frequencies are used for demodulating adding 
and getting the result. This important point of the addition of errors, the fact 
that the spectrum does not have well-defined peaks and it is distributed fairly 
smoothly makes it unwise to use the method of demodulating and adding the sig- 
nals to get a good envelope. 

The Log Approach for Detection of Envelopt^ of Heart Sounds 

Ear is sensitive to log of intensity of the sounds, therefore, in order to get an 
envelope that looks like what the sound sounds like, we can initiate this biologi- 
cal process by electronics. 


Vi 

o 


1 nn 

LOG Vj 





2 LOG V: 



’o 

-►O 


Vj is the rectified (full wave) of the heart sound signal. The Op Amp puts out 
2 log Vj which is log Vj‘ which is log of intensity (pi’oportional to it) and the fil- 
ter smooths the output signal. A lime constant or delay could be introduced at 
the output to simulate the refractory periods of audio sensations of the ear. 

The circuit to get the log of the full wave rectified signal is based upon the prop- 
erty of exponential character of current voltage relationship in a transistor. 




U=a„l„ 

Ics is the reverse saturation 
current of the emitter. 
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Ic = a„ l„ [e f'' 


I = 


Combining equations (a) and (b) we get 


|ia = anic, [ef''-""] 
Rs 


V = M in 

f Rsa„I, 


Converting from logp to logio gives: 


Vo = 2.3 -f- loge 

f ies 


Vo = 2.3 — log.oEin +2.3 — log,o — r 

S S i's “n *es 


atT = 27“C -p = 0.026 V Vo = 0. 06 logio Ejn +K 

Empirically it was found that Vo and Ej have the following relationship: 

Vo = 0.062 logio (Ein ) +0.450 


Vo and Ej are in volts. 
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We tried tMs method, but did not succeed to have a meaningful output. This 
was probably because we should have had extra circuits for frequency compen- 
sation of the operational amplifier. 

The correct circuit and information is given in Application Note AN-261 of 
Motorola Semiconductor Products, Inc. 
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Task m 

PHONOCARDIOGRAPHIC PREPROCESSOR 
COMMENTS 


Methodology 

No formal outline was set up for the project. However, the material is well 
ordered in narrative form in the laboratory notebook. Previous work was re- 
viewed and discussion was frequent with the Department of Clinical Engineering 
personnel. The various methods of intensity detection, such as square law de- 
tection by diodes, temperature devices and frequency shifting, were reviewed 
and discussed in the laboratory notebook. 

Electronic circuits were breadboarded in NASA laboratories and parts were 
supplied by NASA. Various parameters of the circuit were measured and the 
data evaluated. 

Results 


A breadboard circuit consisting of a full wave rectifier with simple resistance 
capacity filtering was demonstrated. This circuit did produce an envelope of 
the phonocardiogram. The envelope, however, seems to have excessive time 
constant and was lacking in defining many of the characteristics of the phono- 
cardiogram. Actual phonocardiographic recording on 1/4" magnetic tape was 
used to develop the circuit. It was interesting to note on this project, that the 
solution tended to be a sophisticated engineering approach, while the users 
(Clinical Engineering) saw a simplified empiric solution. This is a common 
theme in the area of medical instrumentation. This does not mean that instru- 
mentation should not be technically sound, but rather that simplified solutions 
of less accuracy are many times acceptable. 

Conclusion 

The conclusion of the project, althou^ not specifically stated, was that full wave 
rectification with an appropriate time constant circuit is feasible for the detec- 
tion of phonocardiograms. 

Future Application/Expansion 

Work has continued on this project at the Department of Clinical Engineering. A 
circuit called a "box-car detector" was added to the original breadboard circuit 
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of the full wave rectifier. This seems to have less of a "time-constant" prob- 
lem. The signal still needs the appropriate filtering. The filter characteristics 
will be developed using a computer program on sampled data input to simulate 
the response with the results reviewed by physicians and engineei s until an ac- 
ceptable envelope is presented for the appropriate disease categories. 
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